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Claudia Eberlein^ and Robert Zietal^ 

'Dept of Physics & Astronomy, University of Sussex, Palmer, Brighton BNl 9QH, England 

(Dated: February 9, 2008) 

We derive the non-retarded energy shift of a neutral atom for two different geometries. For 
an atom close to a cylindrical wire we find an integral representation for the energy shift, give 
asymptotic expressions, and interpolate numerically. For an atom close to a semi-infinite halfplane 
we determine the exact Green's function of the Laplace equation and use it derive the exact energy 
shift for an arbitrary position of the atom. These results can be used to estimate the energy shift 
of an atom close to etched microstructures that protrude from substrates. 
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I. INTRODUCTION 



An important aim of current experimental cold atom 
physics is to learn how to control and manipulate sin- 
gle or few atom in traps or along guides. To this end 
more and more current and recent experiments deal with 
atoms close to microstructures, most importantly wires 
and chips of various kinds (cf. e.g. If these 

microstructures carry strong currents then the resulting 
magnetic fields and possibly additional external fields of- 
ten create the dominant forces on the atoms, which can 
then by used for trapping and manipulation. However, 
since an atom is essentially a fluctuating dipole, polar- 
ization effects in those microstructures lead to forces on 
the atoms even in the absence of currents and external 
fields, and for atoms very close to them these can be sig- 
nificant. Provided there is no direct wave- function over- 
lap, i.e. the atom is at least a few Bohr radii away from 
the microstructure, the only relevant force is then the 
Casiniir-Polder force [H], which is the term commonly 
used for the van der Waals force between a point-like po- 
larizable particle and an extended object, in this case - 
the atom and the microstructure. 

The type of microstructures that can be used for atom 
chips and can be efficiently manufactured often involve a 
ledge protruding from an electroplated and subsequently 
etched substrate [6| . Here we model this type of system 
in two ways: first by a cylindrical wire (which is a good 
model for situations where the reflectivity of the electro- 
plated top layer far exceeds that of the substrate), and 
second by a semi-infinite halfplane (which is an applicable 
model if the reflectivities of the top layer and substrate 
do not differ by much). While in reality the electromag- 
netic reflectivity of a material is of course never perfect, it 
has been shown by earlier research that interaction with 
imperfectly reflecting surfaces leads to a Casimir-Polder 
force that differs only by a minor numerical factor from 
the one for a perfectly reflecting surface 0, [I]- Thus 
we shall consider only perfectly reflecting surfaces here, 
which keeps the difficulty of the mathematics involved to 
a reasonable level. 



II. ENERGY SHIFT AND GREEN'S FUNCTION 

We would like to consider an atom close to a reflect- 
ing surface and work out the energy shift in the atom 
due to the presence of the surface. As is well known, if 
the distance of the atom to the surface is much smaller 
than the wavelength of a typical internal transition, then 
the interaction with the surface is dominated by non- 
retarded electrostatic forces 5,, ^, ^] . At larger distances 
retardation matters, but at the same time the Casimir- 
Polder force is then significantly smaller and thus hard 
to measure experimentally (see e.g. [13, [HI [13 )■ 
this paper we shall concentrate on small distances that 
lie within the non-retarded electrostatic regime, as these 
lie well within the range of the experimentally realizable 
distances of cold atoms from microstructures @ . 

In order to determine this electrostatic energy shift 
one needs to solve the (classical) Poisson equation for 
the electrostatic potential <& 



(1) 



with the boundary condition that $ = on the surface 
of a perfect reflector. The charge density p would be the 
atomic dipole /x of charges ±g separated by a vector D, 



p(r) = hm^ q ^(5(3)(r - (ro + D)) - S^^\r - Tq) 



(2) 



for a dipole located at rg. Since the dipole is made up 
of two point charges, one can find a solution to the Pois- 
son equation ([1]) via the Green's function G(r, r'), which 
satisfies 



-V2G(r,r') = (5(3)(r _!.') 



(3) 



subject to the boundary condition that it vanishes for all 
points r that lie on the perfectly reflecting surface. For 
the purposes of this paper it is advantageous to split the 
Green's function in the following way, 



G(r,r')- 



1 



(4) 



The first term is the Green's function of the Poisson equa- 
tion in unbounded space, and thus the second term is a 



solution of the homogeneous Laplace equation, chosen in 
such a way that the sum satisfies the boundary conditions 
required of G(r, r'). 

The total energy of the charge distribution is given by 

M 



E = 



1 



d^r p(r)$(r) 



with 



$(r) = — / d^r' G(r,r')p(r') 
£o J 



(5) 
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Thus for a dipole with the charge density ^ the total 
energy is 



E — lim 



q 



D^o\87r£o |r + D- (r + D)| Stteo |r - r| 
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1 



87reo|r + D-r| 87reo |r - (r + D)| 
+ 1^ [Gniro + D, ro + D) - (ro + D, ro)] 



-:^[GH(ro,ro + D)-GH(ro,ro)] ^ . (7) 

The first two terms in this expression are the divergent 
self-energies of the two point charges. The third and 
fourth terms are the energy of the dipole, which also di- 
verges in the limit D — > 0. None of these are interesting 
for us because they are the same no matter where the 
dipole is located. The energy shift due to the presence of 
the surface is given by the remaining four terms, which 
all depend just on the homogeneous solution G_f/(r,r'). 
In the limit D — > these four terms can be written as 
derivatives of Gh, and thus we obtain for the energy shift 
of a dipole /x = qD due to the presence of the surface 



AE^ ^(;,.V)(/x.V') Gff(r,r') 
2eo 



(8) 



r=ro,r'=ro 



where tq is the location of the dipole. 

When applying this to an atom one also needs to take 
into account that for an atom without permanent dipole 
moment the quantum-mechanical expectation value of a 
product of two components of the dipole moment is di- 
agonal. 



(9) 



in any orthogonal coordinate system. This implies that 
for an atom the energy shift due to the presence of the 
surface reads 
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(10) 




FIG. 1: An illustration of the geometry of the dipole near a 
wire. The radius of the wire is R, and the distance of the 
dipole from the centre of the wire is p. 



out the electrostatic Green's function of the boundary- 
value problem for the geometry of this surface. In the 
next two sections we are going to do this for two differ- 
ent surfaces: for a cylindrical wire of radius R and for 
a semi-infinite halfplane. In either case we assume that 
the surfaces are perfectly refiecting, which enforces the 
electrostatic potential $ to vanish there. 



III. NON-RETARDED ENERGY SHIFT NEAR 
A WIRE 



To calculate the energy shift, we first determine the 
Green's function of the Poisson equation in the presence 
of a perfectly reflecting cylinder of radius R and infinite 
length. A standard method of calculating Green's func- 
tions is via the eigenfunctions of the differential operator. 
In order to find a solution to Eq. ^ we solve the eigen- 
value problem 



(11) 



The eigenfunctions V'n(r) must satisfy the same boundary 
conditions as required of the Green's function. Since — 
is Hermitean, the set of all its normalized eigenfunctions 
must be complete. 



Y,\Mr)){Mr')\=S'-'Hr-r') 

n 

Thus one can write 



(12) 



r=ro,r'=ro 



Thus the central task in working out the non-retarded en- 
ergy shift of an atom in the vicinity of a surface is to work 



If we apply this method to unbounded space we can eas- 
ily derive a representation of the Green's function in un- 
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bounded space in cylindrical coordinates, 

1 I ^ , r k 



- y 



4TT\r - r'l 4^2 ^Z^^ ^a. dk ^ 
Performing the k integration [3, 6.541(1.)], we find 



47r|r 
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t = — V 

-r'l 27r2 ^ 



xI„,{kp)K„i{kp') , for p < p' . (15) 



To derive the Green's function (|3|) that vanishes on the 
surface of the cylinder p = R we now just need to find 
the appropriate homogeneous solution G'ij(r, r'), which 
satisfies 



V2GH(r,r') =0 



(16) 



The general solution of the homogeneous Laplace equa- 
tion in cylindrical coordinates can be written as 



E 



m— — oo 



6"""^+"*^ [A{m, K)Im{np) + B{m, K)Km{np)] 

(17) 

where A and B are some constants. Since G(r,r') and 
therefore GHi^, r') must be regular at infinite p we must 
have A{m, k) = 0. This and the requirement that the 
sum of Eq. and Gh{i^, r') must vanish at p = i? lead 
to 



1 ^ fOO 

rn— — OO ^ 
/m(Ki?) 



-,im(0— )+iK(2:— js') 



Km{Kp)K,n{np') . (18) 



The energy shift can now be determined by applying for- 
mula (jlOp in cylindrical coordinates. Taking into account 
the symmetry properties of the modified Bessel functions 
we find for the energy shift of an atom whose location is 
given through (p, (j), z) 



AE : 



--—[E,{pI) +E^{pI) + E^{pI)] (19) 



with the abbreviations 
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7rp2 ^ 
' ^ — 1 
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TO" / d/i . „s 



[K,n{Kp)Y 



an n 1 / -M 



Km(nR) 



[Km(Kp)Y 



The prime on the sums indicates that the m = term is 
weighted by an additional factor 1/2. 
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FIG. 2: The energy shift due to the z component of the atomic 
dipole, multiplied by d^, where d is the distance of the dipole 
to the surface of the wire. The solid line is the exact expres- 
sion calculated numerically, the dot-dashed line is the m = 
term alone, and the dashed line is the m = Q term plus the 
single integral derived through the uniform asymptotic ap- 
proximation for the Bessel functions, Eq. (|20cp . The cross on 
the vertical axis gives the exact value for d — > 0. 
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FIG. 3: Same as Fig. [2] but for the energy shift due to the p 
component of the atomic dipole. The solid line is the exact 
expression calculated numerically, the dot-dashed line is the 
m = term alone, and the dashed line is the m = term plus 
the single integral derived through the uniform asymptotic 
approximation for the Bessel functions, Eq. (|20ap . 



Using numerical integration packages like those built 
into Mathematica or Maple, one can evaluate these con- 
tributions to the energy shift. We show the numerical 
results in Figs. [IHH We have chosen to show the various 
contributions to Eq. as a function of the ratio of the 
distance d — p ~ R oi the atom from the surface of the 
wire to the wire radius R and multiplied by (fi so as to 
plot dimensionless quantities. 

For most values oi d = p — R the integrals over k, con- 
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with the abbreviation 



(20c) 



FIG. 4: Same as Fig. [2] but for the energy shift due to the 
component of the atomic dipole.The solid line is the exact 
expression calculated numerically, the dot-dashed line is the 
m = 1 term alone, and the dashed line is the single integral 
derived through the uniform asymptotic approximation for 
the Bessel functions, Eq. (|20b|) . 



verge quite well, as for large k the dominant behaviour of 
integrands is as exp(— 2Kd). Likewise is the convergence 
of the sums over m very good for reasonably large values 
of d. In fact, as shown by dot-dashed lines in Figs. [2H11 
convergence is so good that from d/i? « 20 upwards it 
is fully sufficient to take just the first summand in each 
sum for S. However, convergence is less good for small d 
and thus the numerical evaluation of the energy shift gets 
more and more cumbersome the closer the atom is to the 
surface of the wire. Neither the integrals over k nor the 
sums over m converge very well, so that it is worthwhile 
finding a suitable approximation for small d. Another 
motivation for a detailed analysis of the limit d — > is of 
course also to check consistency: if the atom is very close 
to the surface of the wire {d <C i?) then the curvature of 
the wire cannot have any impact on the shift any longer 
and the energy shift should simply be that of an atom 
close to a plane surface, which is well known p^ . 



To find a suitable approximation for small d, we sepa- 
rate off the m = terms in 5p and S^. In all the other 
summands and in we scale by making a change of 
variables in the integrals to a new integration variable 
X = np/m. The dominant contributions to those inte- 
grals and the sums come from large x and large m, so 
that one can approximate the Bessel functions by their 
uniform asymptotic expansion [l^, 9.7.7-10]. Then the 
Bessel functions then the sums over m become geomet- 
ric series and can be summed analytically. Taking just 
the leading term in the uniform asymptotic expansions 
for the Bessel functions we find the following approxima- 



A(x) 



1 + Jl + X 



2 



These are easy to evaluate numerically. We show the nu- 
merical values as dashed lines in Figs. [SHU Furthermore, 
these approximations allow us to take the limit d ^ 0. 
Taking this limit under the x integrals and retaining only 
the leading terms in each case, we can carry out the x 
integrations analytically and obtain 



1 

M3 



1 



1 



16d3 



16d3 



(21) 



Inserting these values into Eq. (fT9|) we see that they give 
the energy shift of an atom in front a perfectly reflecting 
plane [la ] . This is an important consistency check for our 
calculation, as the atom should not feel the curvature of 
the surface at very close range. In Figs. [2Hl]the limiting 
values (I^Tj) are marked as crosses on the vertical axes. 
Since we have taken along only first term of the uniform 
asymptotic expansion for each of the Bessel functions, 
we cannot expect any agreement of the approximations 
(HO^c) beyond leading order. However, in practice these 
approximations work quite well beyond leading order: as 
Figs. [2H11 show, approximations (|20bp and (|20cp work 
very well for almost the entire range of d; approximation 
(|20ap works reasonably well for small d and then again 
for large d (because at large distances the m = term 
dominates everything else). In this context we should 
also point out that the m = contributions to Sp and 
Ez do not contribute to leading order, but behave as d~^ 
in the limit d ^ and thus need to be taken into account 
for asymptotic analysis in this limit beyond leading order. 

At large distances d the energy shift (fTO)) is dominated 
by the first terms in each of the sums over m, i.e. the 
m = terms for Sp and and the m = 1 term for 
S^. For large d Sp and behave as l/(d'^lnd), but 
there is no point in giving an asymptotic expression as 
further corrections are smaller only by additional powers 
of 1/ In d and these series converge far too slowly to be 
of any practical use. falls off faster: to leading order 
we get « 37ri?V(32d^) for large d. 
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of summation over eigenfunctions of the operator — V^, 
as explained at the beginning of Section IIIII If we de- 
termine the set of eigenfunctions ipn{j^) that satisfy the 
correct boundary conditions then the Green's function 
([T3| wiU also satisfy these boundary conditions. Normal- 
ized eigenfunctions that vanish at = and (j) — 2tt 
and that are regular at the origin and at infinity are in 
cylindrical coordinates 



z) = 



7= e'"" J™/2 (fcp) ^ sm ^ , (22) 

with the corresponding eigenvalue A„ = + k'^. Con- 
struction (fTS]) then gives the Green's function 



FIG. 5: An illustration of the geometry of a dipole near a 
semi-infinite halfplane. The distance of the dipole from the 
edge is p, and from the surface of the plane it is psin 0. 

Furthermore, if d is not just large compared to R but 
also compared to the typical wavelength of an internal 
transition in the atom, then the energy shift would be 
dominated by retardation effects, which cannot be calcu- 
lated with the electrostatic approach of this paper, but 
which require a quantization of the electromagnetic field 
0. For this reason the extreme large- distance limit of 
the energy shift (fT9)) is not of practical importance and 
thus we see no reason to report on it in more detail. 

IV. NON-RETARDED ENERGY SHIFT NEAR 
A SEMI-INFINITE HALFPLANE 

Next we wish to calculate the energy shift of an atom 
in the vicinity of a perfectly reflecting halfplane. The 
geometry is sketched in Fig. [5l The required boundary 
conditions are that the electrostatic potential $(r) van- 
ishes at the angles = and </> = 27r. We shall calculate 
the Green's function for this case by using the method 



G(r, r') = -^y dk dK ) 

X j!!i{kp)Jm.{kp') sin ^ sin ^ . 

This result is in agreement with the limiting case of the 
Green's function for the electrostatic potential at a per- 
fectly reflecting wedge TT] if the wedge is to be taken to 
subtend a zero angle and extent to an infinite radius. 

As the energy shift (fTU]) depends only on the homoge- 
neous part Gi/(r,r') of the solution, we need, according 
to Eq. (IH), to subtract the free-space Green's function, 
which we have already written down in Eq. (|14p . Noting 
that the free-space Green's function is of course symmet- 
ric under the exchange of and 0', we can write it in the 
form 

1 1 oo rOQ poo 1 

4^|r-r'| 2^2 ^2 + ^2 

X COSm(0 - (/)') Jmikp)J.m{Kp) . 

Taking the difference between G(r, r') above and this ex- 
pression and carrying out the integration over k by clos- 
ing the contour in the complex plane and determining 
the residue, we obtain 



1 \ ^ 

Gff(r,r')--— / dfc e-'^l^-^'l <^ ^ ' J„(M^m(V) [cosm(0 - 0') + cos 

oo 

+ 51 J,n+^ikp)Jm+^ikp' 



m{(j) + (t)')] 



m=0 



COS I m + - 



COS I m + — 



r 



As before, primes on sums over m indicate that the m = dard formulae [ll, 9.1.79], 
term is weighted by an additional factor 1/2. 



The sum in over the product of Bessel functions 
with integer indices can be carried out by applying stan- 



' Jm(kp)J„i{kp') cosm{(j> ± (j)') 



2 MkVp^ + P'^ - ^PP' cos(0 ± </.')) 



(23) 
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Then the k integration can also be carried out in these 
terms by applying well known formulae 11.4.39]. 
Thus we obtain for the terms involving integer indices 
of the Bessel functions 

/•OO ^ 

/ dk e-'^l^-^'l J2 ' Jm,{kp)Jm{kp')cosm{(t)±(t)') 

Jo ^-n 



2 ^[z - Z'Y + p2 + p'2 _ 2pp' cos{(j} ± (j)') 

The terms in with Bessel functions of half-integer 
indices are considerably more difficult to deal with. The 
only relevant formula we could find anywhere is [isl . 
5.7.17.(11.)] 



^ J„+i(fcp)>^m+i(V)cos 



m + - \ a 



sin t 



•y/t^ — fc^p^ — fc^p'^ + 2fc2pp' cos a 



(24) 



with the integration limits 



ti = y/k^p'2 + k^p'^ - 2fc2pp'cosa , t2=k{p + p'). 

Ref. [13] does not give any references, so that the origin 
of this formula cannot be traced. Although the ranges of 
applicability are usually given for formulae in [l^ . they 
are absent in this particular case. However, inspection 



reveals that the formula cannot be valid for the whole 
range < a < 27r but must be restricted to the range 
< a < tt: the right-hand side has a periodicity of 27r, 
i.e. it is the same for a = and a = 27r, but the left- 
hand side differs in sign for these two values of a, as 
cosO = 1 but cos(2to + 1)tt = —1. Thus for the range 
< 10 ± 0'] < TT we apply ^ with a = ± 0', and for 
the range tt < ± (j>' \ < 27r we set a = 2tt ~ {(pzL (/>'). If 
we scale the integration variable from t to s ~ t/k then 
the subsequent k integration is elementary, 



dA: e '"'^ ^'1 sinfcs = 



s2 + (z - Z') 



A2 



(25) 



The remaining integration over s can be carried out by 
changing variables from s to t; = s"^ — sf, so that |14l . 
2.211] 



ds 



{z-z'Y 



2 _ „2 



y/sj + {Z- Z'f 



: arctan - 



- s2 
*2 *1 



sl + iz- z'f 



Along these lines and distinguishing carefully between 
the cases < 10±0'] < tt and tt < l(/)±0'] < 27r we obtain 
the following exact expression for the homogeneous part 
of the Green's function 



Gh{v,v') = -- . ... . , = + 



47r 1 2^(z - z'Y + p2 + p'2 _ 2pp' cos((/> - </)') 2^(z - z'Y + p2 + p'2 _ 2pp' 



cos 



sgn(sin 10 + (/)']) ^ / 2pp' [I + cos{(f) + (f)')] 

arctan - ' 



TTyJiz - z'Y + p2 + p'2 _ 2pp' COs(0 + (j)') \ {z ~ z'Y + P^ + p'^ - 2pp' COs(0 + (j)') 

sgn(sinl0-01) ^^^^^^ / 2pp'[l + cos(0-0O] ~ ] 

7r^(z - z' f + p2 + p'2 _ 2pp' cos((/> - 0') V (■^"^')'+P' + P"-2pp'cos(0-(/)') J 



Calculating the energy shift is now straightforward. 
We apply Eq. (jTU]) and find the exact energy shift of an 
atom located at (p, 0, z) 



Sp(p2)+s^(p2)+S,(p2)] 



(27) 



with the abbreviations 



487rp3 
1 



487rp3 
1 



167rp'^ sin 

^ COS(/) 

Sttp'^ sin^ 
coso!' 



247rp3 I67rp3 sin^ 



(^-0)(1 



IGttp'^ sin'^ ■ 



167rp'^ sin'^ ( 



167rp-^ sin 



Here the applicable range of is < < tt. Geometries 
with (p in the range tt < < 27r are obviously just mirror- 
images of those with < < tt, so that one could simply 
replace by 27r — 0. 



7 



An important cross-check is the hmit cj) ^ 0. Very 
close to the surface the edge of the half-plane should not 
affect the energy shift, as the distance of the atom to the 
edge is p but the distance to the surface is d « pt/i <C p. 
Thus the energy shift should be the same as that in front 
of a plane. Indeed, the leading terms in the limit (j) —^ 
are 

which, if inserted into ([?/)) ■ give the energy shift of an 
atom in front of an infinitely extended reflective plane a 
distance p(j) away. Note that in contrast to Eq. (pij) for an 
atom very close to a cylindrical surface, the component 
of the dipole that is normal to the surface is now /i^. 

V. SUMMARY 

We have calculated the non-retarded Casimir-Polder 
force on a neutral atom that is either a distance d away 
from a cylindrical wire of radius R, or somewhere close 
to a semi-infinite halfplane. In both cases we have, for 
simplicity, restricted ourselves to calculating the force for 
the atom interacting with a perfectly reflecting surface. 
We have worked in cylindrical coordinates (p, (j), z) and 
chosen the z axis along the centre of the wire and along 
the edge of the halfplane, respectively. The energy shifts 
in a neutral atom due to the presence of the reflecting 



surface nearby depend on the mean-square expectation 
values of the dipole moments of the atom along the three 
orthogonal directions. For an atom close to a reflecting 
wire the shift is given by Eq. (jl9p . and for an atom near 
a semi- infinite halfplane by Eq. ([27)1 . The problem of the 
atom and the wire has two independent length scales, 
the radius of the wire and the distance of the atom from 
the wire, and accordingly the shift varies with the ratio 
between them. We have provided analytical and numer- 
ical approximations for both small and large values of 
this ratio in Section III. By contrast, in the problem of 
an atom close to a semi-infinite halfplane the distance of 
the atom to the surface is the only available length scale, 
so that the energy shift does not depend on any further 
parameters. This suggests that one could find an exact 
expression for the energy shift, as we have indeed man- 
aged to do. It is nevertheless surprising that this exact 
expression, Eq. (j27p . can be given in terms of elementary 
functions and is so simple. Both model situations can be 
useful for estimating the energy shift of an atom close to 
microstructures that consist of a ledge and possibly an 
electroplated top layer. 
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